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correspond to an integrated luminosity of 5.0 fb^^ of pp collisions at a centre-of-mass 
energy of 7 TeV collected with the CMS detector at the Large Hadron Collider. No 
excess of events is observed beyond the expected standard model background and 
therefore upper limits are set on the square of the mixing element, | V^nP/ for £ = e,ji, 
as a function of Majorana neutrino mass. These are the first direct upper limits on the 
heavy Majorana-neutrino mixing for > 90 GeV. 
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1 Introduction 

The non-zero masses of neutrinos, confirmed from studies of their oscillations among three 
species, provide the first evidence for physics beyond the standard model (SM) [H. The small- 
ness of neutrino masses underscore the lack of a coherent formulation for the generation of 
mass of elementary particles. The leading theoretical candidate for accommodating neutrino 
masses is the so-called "seesaw" mechanism 121-0, where the smallness of the observed neu- 
trino masses (niy) is attributed to the largeness of a mass (otn) of a new massive neutrino state 
N, with niv ~ yl'i^^/ ^N, where i/y is a Yukawa coupling of v to the Higgs field, and v is the 
Higgs vacuum expectation value in the SM. Owing to the new heavy neutrino's Majorana na- 
ture, it is its own antiparticle, which allows processes that violate lepton-number conservation 
by two units. Consequently, searches for Majorana neutrinos are of fundamental interest. 

The phenomenology of searches for heavy Majorana neutrinos at hadron colliders has been 
considered by many authors II6 HT3I . Our search follows recent studies at the Large Hadron 
Collider (LHC) [TV. ^T2] that use a model-independent phenomenological approach, with 
and as free parameters, where V^n is a mixing element describing the mixing between 
the heavy Majorana neutrino and the SM neutrino of flavour i. Previous direct searches 
for heavy Majorana neutrinos based on this model have been reported by the L3 |!T4J and 
DELPHI |il5;l Collaborations at the Large Electron-Positron Collider. They have searched for 
Z — >• VfN decays and set limits on jV^Np ss a function of jwn for Majorana neutrino masses 
up to approximately 90 GeV. The ATLAS Collaboration at the LHC has also reported limits on 
Majorana neutrino production [ l6l HZl in the context of an effective Lagrangian approach ITSl 
and the Left-Right Symmetric Model ||T9ll20| . Indirect limits on | VeN 1^ have been obtained from 
the non-observation of neutrinoless double beta decay [21], resulting in 90% confidence level 
(CL) limits of | VgNp/ wJn < 7 x 10^^ TeV^^. Precision electroweak measurements can also be 
used to constrain the mixing elements resulting in indirect 90% CL limits of [FeNp < 0.0066, 
|y,,Np < 0.0060, and [y^Np < 0.016 [22J. 

We report on a search for the production of a heavy Majorana neutrino in proton-proton col- 
lisions at a centre-of-mass energy oi ^ys = 7 TeV at the LHC using a set of data of integrated 
luminosity 4.98 ± 0.11 fb^^ collected with the Compact Muon Solenoid (CMS) detector. The 
principal Feynman diagram for this process is shown in Fig. [l] The Majorana neutrino can de- 
cay to a lepton with positive or negative charge, leading to events containing two leptons with 
the same or opposite sign. Same-sign events have much lower backgrounds from SM processes 
and therefore provide an accessible signature of Majorana neutrino production. We search for 
events with two isolated leptons of same sign and flavour or e^e^) and at least two ac- 

companying jets. Contributions from SM processes to such dilepton final states are very small 
and the background is dominated by processes such as multijet production, in which leptons 
from b-quark decays or from jets are misidentified as isolated prompt leptons. 

2 Signal simulation and data selection 

The production and decay process is simulated using the event generator described in Ref. 
and implemented in ALPGEN ||23|| . We use the CTEQ5M parton distribution functions Il24ll . Par- 
ton showering and hadronization are simulated using PYTHIA II25I . The Monte Carlo generated 
events are interfaced with CMS software, where Geant4 [1261 detector simulation, digitization 
of simulated electronic signals, and event reconstruction are performed. Monte Carlo sim- 
ulated events are mixed with multiple minimum bias events with weights chosen using the 
distribution of the number of reconstructed primary vertices observed in data to ensure correct 
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2 Signal simulation and data selection 




Figure 1: The lowest-order Feynman diagram for production of a Majorana neutrino N. The 
charge-conjugate diagram also contributes and results in a £^£^qq' final state. 

simulation of the number of interactions per bunch crossing (pileup). The cross section for the 
process shown in Fig. [I] for | Vij^ P = 1 has a value of 866 pb for = 50 GeV, which drops to 
2.8 pb for = 100 GeV, and to 83 fb for wn = 210 GeV Hill . 

The CMS detector is described in detail in Ref. [27]. Its central feature is a superconducting 
solenoid, which provides a magnetic field of 3.8 T along the direction of the counterclockwise 
rotating beam, taken as the z axis of the detector coordinate system, with the centre of the 
detector defined to be z = 0. The azimuthal angle (p is measured in the plane perpendicular to 
the z axis, while the polar angle 6 is measured with respect to this axis. Muons are measured 
in four layers of gaseous ionization detectors embedded in the steel return yoke of the magnet, 
while all other particle detection systems are located inside the bore of the solenoid. Charged 
particle trajectories are measured in a silicon pixel and strip tracker covering < ^ < 27r in 
azimuth and |f/| < 2.5, where rj is the pseudorapidity defined as = — ln[tan(0/2)]. The 
tracker is surrounded by a lead tungstate crystal electromagnetic calorimeter (ECAL) and a 
brass/ scintillator hadron calorimeter that are used to measure the energy of electrons, photons, 
and hadronic jets. A two-level trigger system selects the most interesting events for analysis. 

Dilepton triggers are used to select the signal sample. Depending on the average instantaneous 
luminosity of the LHC, dimuon events are recorded using triggers requiring the presence of 
two muons with transverse momenta (px) either above 7 GeV for both muons, above 13 GeV 
for one muon and above 8 GeV for the second, or above 17 GeV for one muon and above 8 GeV 
for the second. Dielectron events are recorded using triggers that require one electron with 
transverse momentum above 17 GeV, and above 8 GeV for the second. Trigger efficiencies are 
measured using Z — > and Z — > e+e^ events selected in data, and are found to be (96.0 ± 

2.0)% for muons and (98.5 ± 1.0)% for electrons. 

Additional selections are performed offline to ensure the presence of well-identified muons, 
electrons, and jets. Events are first required to have a well-reconstructed primary vertex based 
on charged tracks reconstructed in the tracking detectors. 

Muon and electron candidates are required to have \t]\ < 2.4 and to be consistent with orig- 
inating from the primary interaction vertex. Muon candidates are reconstructed by matching 
tracks in the silicon tracker to hits in the outer muon system, and are also required to sat- 
isfy specific track-quality and calorimeter-deposition requirements. Electron candidates are 
reconstructed from energy depositions in the ECAL. These are matched to tracks in the sili- 
con tracker and are required to satisfy shower distribution and cluster-track matching criteria. 
Electron candidates within AR = ^ (A?/)^ + (A^)^ < 0.4 of a muon candidate are rejected to 
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remove spurious electron candidates formed from the track of a muon that has an associated 
photon from bremsstrahlung. Electron candidates from photon conversions are suppressed by 
looking for a partner track and requiring that this track has no missing hits in the inner layers 
of the silicon tracker. 

Electron and muon candidates must be isolated from other activity in the event by requiring 
their relative isolation (Zrei) to be less than 0.1. Here Jj-ei is defined as the scalar sum of transverse 
track momenta and transverse calorimeter energy depositions present within AR < 0.3 of the 
candidate's direction, excluding the candidate itself, divided by its transverse momentum. 

Jets and the missing transverse energy in the event are reconstructed using the objects defined 
in the particle-flow method 1281 129II . Jets are formed from clusters based on the anti-fcx algo- 
rithm [30 1, with a distance parameter of 0.5 and are required to be within the pseudorapidity 
range \fj\ < 2.5 and to have transverse momentum pj > 30 GeV. The missing transverse energy 
is required to be less than 50 GeV. 

Events in the muon channel are required to contain two same-sign muons, one with pj greater 
than 20 GeV and the other with pj greater than 10 GeV. Events with an opposite-sign third 
muon that combines with one of the other candidate muons to give a invariant mass 

within the window for a Z boson of 76-106 GeV are excluded. In the electron channel, events 
are required to contain two same-sign electrons, one with pi greater than 20 GeV and one with 
Pj greater than 10 GeV. Events containing any third electron candidate are rejected. Overall sig- 
nal acceptance includes trigger efficiency, geometrical acceptance, and all selection criteria. In 
the muon channel, the overall acceptance for heavy Majorana neutrino events ranges between 
0.43% for = 50 GeV to 29% for = 210 GeV. For the electron channel, the corresponding 
efficiency changes from 0.40% to 21% for these masses. The lower acceptance at low is due 
to the smaller average pi of the jets and leptons in these events. 

3 Background estimation 

There are three potential sources of same-sign dilepton backgrounds. The first and most im- 
portant originates from events containing leptons from b-quark decays or generic jets that are 
misidentified as leptons. Examples of this background include: (i) multijet production in which 
two jets are misidentified as leptons; (ii) W(— ?► Iv) + jets events in which one of the jets is 
misidentified as a lepton; and (iii) tt decays in which one of the top quarks decays giving a 
prompt isolated lepton (t Wb (.Vfb), and the other lepton of same charge arises from a b- 
quark decay. From Monte Carlo studies we find the dominant contribution to this background 
is from multijet production, with the sum of W(— >• Iv) + jets and tf events comprising approxi- 
mately 15-35% of the total misidentified lepton background. These backgrounds are estimated 
using control samples in collision data as described below. 

To estimate the misidentified lepton background, an independent data sample enriched in mul- 
tijet events is used to calculate the probability for a jet that passes minimal lepton selection 
requirements to also pass the more stringent requirements used to define selected leptons. The 
lepton candidates passing the less stringent requirements are referred to as "loose leptons" 
and their misidentification probability is calculated as a function of transverse momentum and 
pseudorapidity. This probability is used as a weight in the calculation of the background in 
events that pass all the signal selections except that one or both leptons fail the tight criteria. 
This sample is referred to as the "orthogonal" sample. 

The misidentification probability is applied to the orthogonal sample by counting the number 
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3 Background estimation 



of events in which one lepton passes the tight criteria, while the other lepton fails the tight se- 
lection but passes the loose selection (Nnn), and the number of events in which both leptons fail 
the tight selection, but pass the loose criteria (Nnn)- The total contribution to the signal sample 
(i.e. the number of events when both leptons pass the tight selection, Nnn), is then obtained by 
weighting events of type nn and tm by the appropriate misidentification probability factors. To 
account for double counting we correct for Wn events that can also be nn. 

In the muon channel, loose muons are defined by relaxing the muon isolation requirement from 
Jj-ei < 0.1 (used to select signal events) to Jj-ei < 0.8. In the electron channel, loose electrons are 
defined by relaxing the isolation from 0.1 to 0.6, and by removing a requirement on transverse 
impact parameter normally used for tight electrons. 

We evaluate the method used to estimate the background from misidentified leptons by check- 
ing the procedure using Monte Carlo simulated event samples in which the true background 
is known. The misidentification probabilities are obtained from multijet events and are used to 
estimate the misidentified lepton backgrounds in tt, W + jets, and multijet events by applying 
the background estimation method described above. The differences between the estimated 
backgrounds and the true number of events in the Monte Carlo samples is used as input to the 
overall systematic uncertainty. 

The overall systematic uncertainty on the misidentified lepton background is determined from 
the variation of the background estimate with the loose lepton definition and the variation 
with the leading jet pj requirement in the data samples used to measure the misidentification 
probability, as well as from the independent Monte Carlo validation studies described above. 
The overall uncertainty is found to be 35%. 

The second contribution to the background is from £^£^ events in which the charge of one of 
the leptons is mismeasured. Since the charge mismeasurement probability is very small for 
muons, this background is significant only in the electron channel. The background from mis- 
measurement of electron charge is estimated through probabilities calculated in Monte Carlo 
studies of e+e^ events that pass all selections for signal, except the same-sign requirement. The 
average electron mismeasurement probability is found to be (3.3 ± 0.2) x 10^^ in the ECAL 
barrel region (|j/| < 1.5) and (2.9 ± 0.1) x lO^^ in the ECAL endcap region (1.5 < \rj\ < 2.5). 
To validate the charge mismeasurement probability, we select a control sample of Z — > e^e^ 
events in data, requiring an e^e^ invariant mass between 76 and 106 GeV. We use the differ- 
ence between the predicted and the observed number of e^e^ events, including uncertainties, 
to set a systematic uncertainty of 25% on the background from charge mismeasurement. 

The third background source is the irreducible background from SM production of two gen- 
uine isolated leptons of the same sign, which can originate from sources such as ZZ, WZ, and 
W7 diboson production, ttW, double W-strahlung W^W^qq processes, and double-parton 
scattering (two qq' — > W). These have relatively small cross sections, and are consequently 
estimated using Monte Carlo simulations. We use pythia to simulate ZZ and WZ production 
and MadGraph Il3l1 for the remaining processes. 

After applying all the final selections we observe 65 events in data in the muon channel and 
expect a total SM background of 70 ± 4 (stat.) ± 22 (syst.) events, with the dominant contribu- 
tion of 63 ± 4 (stat.) ± 22 (syst.) events arising from the misidentified muon background. The 
data are in agreement with the estimated background. In the electron channel we observe 201 
events in data, and estimate the total SM background as 219 ± 6 (stat.) ± 62 (syst.) events, with 
the dominant contribution of 177 ± 5 (stat.) ± 62 (syst.) events arising from the misidentified 
electron background. The data are again in agreement with the estimated background. The 
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final estimates for the two channels are given in Table [T] 

Table 1: Observed event yields and estimated backgrounds in the muon and electron channel. 
The first quoted uncertainty on each background estimate is statistical and the second is sys- 
tematic. For the irreducible SM backgrounds, the first uncertainty is due to the statistical error 
associated with the finite size of the Monte Carlo event samples used. 



Source 


^^^^ 




e±e± 


Irreducible SM backgrounds: 








WZ 


3.2 ±0.3 ±0.2 




4.9 ±0.3 ±0.3 


ZZ 


1.0 ±0.1 ±0.1 




2.1 ±0.1 ±0.1 


W7 


0.75 ± 0.27 ±0, 


.07 


1.7 ±0.4 ±0.2 


ttw 


1.06 ± 0.05 ±0, 


.53 


0.62 ± 0.04 ± 0.31 


W+W+qq 


0.76 ± 0.06 ± 0, 


.38 


0.73 ± 0.07 ± 0.37 


W-W-qq 


0.45 ± 0.03 ± 0, 


.23 


0.27 ±0.02 ±0.13 


Double-parton W+W+ 


0.07 ± 0.02 ±0, 


.04 


0.19 ±0.03 ±0.10 


Total irreducible SM background 


7.3 ± 0.4 ± 0.7 




10.6 ±0.6 ±0.6 


Charge mismeasurement background 


Q+0.2 




31.9 ±2.7 ±8.0 


Misidentified lepton background 


63.1 ±4.2 ±22, 


.1 


176.8 ±4.7 ±61.9 


Total background 


70 ± 4 ± 22 




219 ± 6 ± 62 


Data 


65 




201 



4 Systematic uncertainties 

The sources of systematic uncertainty associated with signal efficiency and background esti- 
mates can be summarized as follows. 

(1) The systematic uncertainty on the integrated luminosity is 2.2% ||32| . 

(2) The systematic uncertainty from choice of parton distribution functions is estimated 
from Monte Carlo simulations following the PDF4LHC recommendations [33j and 
is found to be 6% of the signal yield. 

(3) The hard-scattering scale in the alpgen Monte Carlo generator is varied from the 
nominal value of to 4Q^ and Q^/4. The resulting uncertainty on the signal yield 
is 1%. 

(4) The jet energy scale is changed by its estimated uncertainty l|34| resulting in a sys- 
tematic uncertainty of between 3.3% for high Majorana neutrino mass (nif^ = 210 GeV) 
and 14.2% at low mass (m^ = 50 GeV). For low mass events the jets from the Majo- 
rana neutrino decay have lower average pj, leading to larger uncertainties. 

(5) The systematic uncertainty due to the uncertainty on the jet energy resolution 13411 is 
determined to be 0.2% — 1%, depending on m^. 

(6) The uncertainty in modeling event pileup is studied in the Monte Carlo simulations 
and found to be 1%. 

(7) The systematic uncertainty on the estimate of misidentified lepton background is 
35%. 

(8) The systematic uncertainty on the background from mismeasurement of electron 
charge is 25%. 
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6 Summary 



(9) The systematic uncertainties on normalizations of irreducible SM backgrounds are: 
6% for WZ and ZZ 135]; 10% for W7 ISSl; and 50% for the other processes, deter- 
mined by varying the scale and parton distribution functions in Monte Carlo 
simulations. 

In the muon charmel the systematic uncertainty due to the muon trigger, as indicated above, is 
based on studies of Z — > events in collision data, and is determined to be 2% per muon. 

The offline muon selection efficiency is taken from Monte Carlo simulation, and cross-checked 
with data using studies of Z — > ji^ji^ events. The efficiencies measured in data and Monte 
Carlo simulation are found to be in agreement within uncertainties; the systematic uncertainty 
associated with the small differences is 2%. The overall systematic uncertainty due to the muon 
trigger and selection is 4%. In the electron channel the systematic uncertainty from the trigger 
and electron selections is determined in a similar way, and found to be 10%. 

5 Results and discussion 

We see no evidence for a significant excess in the data beyond the backgrounds predicted from 
the SM and set 95% CL exclusion limits on the square of the Majorana-neutrino mixing element 
as a function of m^, using the CLs method [36:-:38J. In the muon channel analysis we set limits 
on |Vj,Np a function of m^, under the assumption [VgNp = |VVnP = 0- In the electron 
channel analysis we set limits on |VeNp as a function of under the assumption |V^nP = 
I = 0- Figure |2|a) shows the resulting upper limits in the muon channel (| Vj,nP vs. /Wn), 

while Fig. |2|b) shows the upper limits in the electron channel (| VeNp vs. m^). These are the 
first direct upper limits on the heavy Majorana-neutrino mixing for > 90 GeV. 

For low /«N the limits in both channels are less stringent than the existing limits from DELPHI 
and L3 shown in Figs. |2]^a) and |2|b), due to the higher backgrounds at the LHC. However, the 
DELPHI and L3 limits are derived from Z — )■ V(N and are restricted to masses below approx- 
imately 90 GeV. The limits reported here extend well beyond this mass. For = 90 GeV we 
find |y^N|^ < 0.07 and l^eNp < 0.22. At = 210GeVwefind |V,,nP < 0.43, while for |FeNp 
the limit reaches 1.0 at a mass of 203 GeV. 

6 Summary 

A search for heavy Majorana neutrinos in and e^e^ events has been performed using a 

set of data corresponding to 5.0 fb^^ of pp collisions at a centre-of-mass energy of 7TeV. No 
excess of events beyond the standard model background prediction is found. Upper limits at 
the 95% CL are set on the square of the Majorana-neutrino mixing element, | V^^P/ for £ = e,}i, 
as a function of Majorana neutrino mass, as shown in Figs.|2]ja) and|2]^b). For = 90 GeV the 
limits are |F,,n|^ < 0.07 and IVeNp < 0.22. At = 210 GeV the limits are | V^^nP < 0.43, while 
for |VeNp the limit reaches 1.0 at a mass of 203 GeV. These are the first direct upper limits on 
the heavy Majorana-neutrino mixing for > 90 GeV. 
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Figure 2: Exclusion region at 95% CL in the square of the Majorana-neutrino mixing element as 
a function of the Majorana neutrino mass: (a) | V^nP vs m^; (b) | VeNp vs. m^. The long-dashed 
black line is the expected upper limit, with one and two standard-deviation bands shown in 
dark green and light yellow, respectively. The solid red line is the observed upper limit, and is 
very close to the expected limit such that the two curves almost overlap. Also shown are the 
upper limits from L3 | [T4| and DELPHI IITSl . The regions above the exclusion lines are ruled out 
at 95% CL. 
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